Towards efficient algorithm deciding separability of distributed quantum states. 
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It is pointed out that separability problem for arbitrary multi-partite states can be fully solved by 
a finite size, elementary recursive algorithm. In the worse case scenario, the underlying numerical 
procedure, may grow doubly exponentially with the state's rank. Nevertheless, we argue that for 
generic states, analysis of concurrence matrices essentially reduces the task of solving separability 
problem inmxn dimensions to solving a set of linear equations in about (j nn +^>- ls j variables, where 
D decreases with mn and for large mn it should not exceed 4. Moreover, the same method is also 
applicable to multipartite states where it is at least equally efficient. 



I. INTRODUCTION 

Presence of quantum correlations in a distributed sys- 
tem is probably the most clear and fundamental marker 
of the system's non-classicality. In applications, quan- 
tum correlations make grounds for the emerging field of 
quantum information technology. In either context, be it 
fundamental or utilitarian, the notion of quantum entan- 
glement plays the key role 0. Not surprisingly then, a 
search for an unambiguous price-tag for the entanglement 
present in a given arbitrary mixed distributed state has 
been one of the important tasks of physics of information. 
It appears, however, that this task is notoriously difficult 
even in the seemingly simple case of bipartite states. A 
yet simpler problem of a finite, operational necessary and 
sufficient separability test for bipartite states in arbitrary 
dimensions is still difficult. 

The claim that an arbitrary bipartite n ® m state g 
is separable means that this state can be expressed as a 
convex combination of projectors on product directions 
in the bipartite Hilbert space TLab ■ This can be formally 
expressed as the following condition: 

The statement 



i = l.. 

u=l 



. 7i, j = 1 ... to, i-j. = l... (mn) 
2-1 = QiM 



1) 



(1) 
is true. 

To decide the truth value of statement Q is usually 
difficult, nevertheless it is always possible in principle. 
The problem to solve is an example of a decision problem 
for the existential theory for reals with s = (mn) 2 — 1 
(number of real parameters in g) quartic (degree d = 4) 
polynomials in k = 2(m + n)(m 2 n 2 — 1) real variables. 
As such, the problem represents a special case of the 
quantifier elimination problem. The fact that there ex- 
ist algorithms solving these problems was first proved by 
A. Tarski |2j. Complexity of his algorithm was, however, 
not elementary recursive. It means that it could not be 
bounded from above by any power of exponentials of a 
finite hight. 

Although considerable progress in the area has been 
made over the years and there are working algorithms 



for the quantifier elimination problem implemented into 
such programs like MAPLE or Mathematica, the exist- 
ing algorithms are far too inefficient to solve non-trivial 
separability problems. The best general purpose algo- 
rithm which we have found is by Basu, Pollack and Roy 
and it solves the decision problem for the existential 
theory using s k + 1 ( f ) (k) arithmetic operations. For sepa- 
rability of a 3 x 3 system formulated as in statement Q 
this means an astronomical number of operations, of the 
order of 81 961 • 4°( 96 °) « 10°( 2412 ) . 

The line of research aiming at a finite operational (in 
principle analytic) necessary and sufficient criterion of 
entanglement in arbitrary distributed states relies on the 
hope that the head on attack on the problem described 
above is unnecessarily complicated. Indeed, even for a 
pair of qubits, solving the decision problem Q directly 
would require 16 121 • 4°( 120) w 10°( 22 °) operations to 
determine entanglement. One knows, however, that in 
this case as well as for qubit-qutrit pairs, a simple par- 
tial transpose criterion solves the problem efficiently Q. 
Moreover, in case of two qubits, Wootters found an effi- 
cient finite algorithm to determine an optimal decompo- 
sition of an arbitrary 2-qubit state [5j . These results sug- 
gest some natural approaches to the separability problem 
in arbitrary dimensions. 

One may thus seek for a method which would reduce 
separability to an eigenvalue problem, like in the partial 
transposition (PT) test [fj. Among the possible exten- 
sions and derivatives of PT, the hierarchy of the PT tests 
as well as the matrix realignment criterion [t| and 
its extensions |n| have the desired, operational (ie. 
eigenvalue-like) form. Moreover, these tests are generally 
more powerful detectors of entanglement than PT. Nev- 
ertheless, they still do not guarantee any universal neces- 
sary and sufficient condition for separability. Entangle- 
ment witnesses method 0, and its optimizations jT^] 
(important form experimental point of view) have similar 
drawbacks, although here, an exhaustive search over all 
the possible pure separable states in principle solves the 
problem. 

To avoid searching over an uncountable set, Gurvits 
reformulated the separability problem as a weak mem- 
bership problem and showed that even in this limited 
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guise, the problem is still NP-hard. In the development 
which followed, Ioannou at al. proved the existence of 
an, in principle, efficient algorithm solving the underly- 
ing weak membership problem |l5j . The structure of the 
algorithm was, however, complicated, which led to tech- 
nical difficulties with its implementation. 

Hulpke and BruB 

showed that irrespective of the 
required accuracy, the search may be limited to the states 
with rational coordinates in a given basis. This result 
combined with the symmetric extensions separability cri- 
terion by Doherty et al. @J guaranteed a finite stopping 
time (although no upper bound on it) for a separabil- 
ity test on all the states but those exactly on the border 
between separable and entangled. 

Wootters' solution to the separability of two qubits re- 
lied on the concept of concurrence and its generalization 
for mixed states, concurrence matrix [{|. In particular, 
introduction of concurrence matrix allowed him to derive 
a simple closed formula for entanglement of formation in 
terms of the singular values of the matrix. After several 
attempts (see the concept of concurrence matrix 

was extended to arbitrary bipartite systems [l^. There, 
it appeared that a bilinear combination of concurrence 
matrices, a biconcurrence matrix maybe the fundamen- 
tal object to investigate in connection with separability 
of bipartite quantum stats. This line of thought was later 
confirmed by Mintert, Kus and Buchleitner who used bi- 
concurrence to formulate powerful tests for entanglement 
of bipartite systems • Recently, the same authors suc- 
cessfully generalized biconcurrence to cover multipartite 
states as well |2(|. This led to powerful tests for differ- 
ent kinds of multi-partite entanglement. However a tight 
necessary and sufficient operational test still remained 
out of reach, since it would usually need optimization 
procedures. 

In what follows, we reexamine the role of biconcurrence 
in solving separability problem for arbitrary distributed 
states. In particular, we study the underlying decision 
problem formulated in terms of the eigenvectors of bi- 
concurrence matrix and show that not only for low rank 
density matrices (like it was in the analysis of range cri- 
terion 0,H3), but also for generic density matrices in 
many dimensions the resulting solutions are computa- 
tionally manageable. 

In exceptional cases, the resulting numerical calcula- 
tions may still grow super exponentially with the increas- 
ing rank(g). Therefore we do not claim a complete so- 
lution to the separability problem. Nevertheless, we can 
quote evidence that for generic multipartite states the 
calculations grow polynomially with the size of the in- 
vestigated system. An improvement of the method which 
would entirely remove the difficult exceptional cases can 
hardly be expected since separability problem appears to 
be computationally hard [f4| . 



II. BIPARTITE SEPARABILITY AND 
CONCURRENCE. 

Entanglement of pure bipartite states is well under- 
stood, so let us begin by considering such a state 

\>I>)ab =^2a i:i \i) A ® \j) B (2) 

id 

(the i's and the j's here number orthonormal basis vec- 
tors in the Hilbert spaces held by Alice (A) and Bob (B) 
respectively). Separability of \iP)ab means that the right 
hand side of (J2J factorizes, i.e., that 

H)ab = ^2fi\i)A ®9j\j)B (3) 

id 

This is equivalent to the following conditions for the ma- 
trix elements aij 

V(i, j, k, I) a tj a u - auctkj = (4) 

For a pure state in m x n dimensions, equation Q im- 
poses N — m "( m ~ 1 )(" 1 ) non-trivial conditions. They 
correspond to the requirement that all the 2x2 minors 
of matrix [a] are zero. The degree to which a given pure 
state violates condition J3J can thus be measured by the 
length of the A/"-dimensional vector C with components 

C a = Ci/\k,jAl — 2(aija,ki — a-aUkj) (5) 
This length, 



coincides with concurrence defined for pure multi- 
dimensional systems by P. Rungta et al. 17]. 

Separability of a mixed state g means that there ex- 
ist a decomposition of g, where all the contributing pure 
states have vanishing concurrence I©- One should no- 
tice here that the decomposition of a given state is not 
unique. In fact, any two ensembles of (sub-normalized) 
pure states { IV'^) } an d {1^)} realize the same state 
g = J^fj, I^XV^I iff they are related by a unitary trans- 
formation |2^] 

\^) = Y J u lw \r) (7) 

V 

To decide existence of a separable decomposition of a 
given state is in general difficult. One of the sources of 
this difficulty may be associated with the fact that the 
set of pure state concurrences contributing to a mixed 
state is not closed with respect to the changes of the de- 
composition. The objects nearest to C which transform 
well with the decomposition changes are concurrence ma- 
trices. For a given state decomposition 

= £KWI (8) 
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with sub-normalized pure components = 
J2i j a ijV) 'A ® one can define N symmetric 

concurrence matrices 

CT = a?y kl - a»a v kj + <•< - a>£. (9) 

When the state decomposition is changed from 
to {1^)} via transformation 10, then the concurrence 
matrices C£ v change into 

CT = J2 U ^CfU^ (10) 

So, in terms of concurrence matrices, separability means 
that there exists such a state decomposition for which 
the diagonals of all the TV concurrence matrices C a are 
zero. One can immediately notice that in this context two 
qubits represent an exceptional case: there is only one 
concurrence matrix. This simplification allowed Woot- 
ters to produce an elegant solution of the separability 
problem for two qubits 5]. Moreover, his solution re- 
sulted in a simple and efficient algorithm for the deter- 
mination of an optimal decomposition of both separable 
and entangled two-qubit states. 

The 2-qubits solution does not, however, generalize 
easily for systems in higher dimensions. There, in terms 
of concurrence matrices, a bipartite state is separable iff 
there exists isometry U such that 

= u^cfu^ 13 = o (ii) 

a/3 

for all values of a simultaneously. 

It appears that the question of existence of U which 
satisfies the whole set of conditions (|11|) is much more 
difficult to answer than that of existence of U which sat- 
isfies only one of these conditions (for a single given a) 
solved in 0. In particular, in a state with bound entan- 
glement, none of the C CT 's will show entanglement sepa- 
rately, nevertheless, regardless of the choice of the local 
bases, it will be impossible to find a decomposition which 
zeros the diagonals of all the C„'s together. 

To make the problem independent of the local choice 
of bases and possibly reduce its complexity, one may con- 
sider a single object, biconcurrence instead of the set of 
concurrence matrices For clarity, one may permute 
indexes in the original definition (formulae (22) and (23) 
in 0|) and define biconcurrence B as 

Bfi v; a f3 = J2c%"C* a P (12) 

(the asterisk denotes complex conjugate). Written in 
this form, B represents a manifestly positive operator in 
the vector space of symmetric concurrence matrices. 

In terms of biconcurrence, a state g is separable iff 
given an initial decomposition of the state, e.g., the eigen- 
decomposition, there exists isometry U such that 

G = U^U^B^ ^(jj ia U ifi )* = (13) 

i— 1 ^iuaj3— 1 



This, together with the conditions making U represent 
an isometry, allows to associate entanglement with strict 
positivity of a real quartic form in at most 2mn{{mri) 2 — 
1) variables. The latter problem was formally solved by 
Jamiolkowski back in 1972 [24|. The level of complexity 
of this solution was, however, to high to be of practical 
relevance in the present context. 

On the other hand, one may notice that existence of 
an isometry U solving (|13J) is equivalent to the existence 
of such U which zeros the diagonals of all those eigen- 
vectors of B (hereafter denoted by Tjf u ), which belong 
to non-zero eigenvalues. By using biconcurrence in this 
way, Mintert, Kus and Buchleitner were able to derive 
simple, powerful lower bounds on the entanglement of 
formation for arbitrary bipartite states |19|. In particu- 
lar, it appeared that unlike the C CT 's, single matrices T CT 
can indicate bound entanglement, e.g., in the entangled 
states introduced in [23] and showed to be non-distillable 
in [H. 

Concurrence matrices can be powerful indicators of en- 
tanglement even when it is not seen in any single linear 
combination of T CT matrices, i.e., when the diagonal of 
every linear combination of T a 's can be brought to zero 
separately. To see it, it is enough to notice that the equa- 
tion set 111(1 (or an equivalent equation set with the T CT 's 
substituted for C CT 's) represents a set of up to JV homoge- 
neous quadratic equations for r elements of the /z'th row 
of isometry U and that the elements u a of each row sat- 
isfy the same set of equations (when one begins with the 
eigendecomposition then r is equal to the state's rank). 



^ t;V«^o, (7 = 12,..,^ (w) 

a,/3=l 

When the solution set of this equation set is empty 
than one can clearly claim entanglement right a way. 
Otherwise, whenever the number of equations exceeds 
r — 2, one may expect a finite number of solutions, i.e., 
a finite number of possible rows of U (modulo normal- 
ization). One should expect this to be the prevailing 
situation since a generic state in m x n dimensions is of 
rank r g = mn while the corresponding number of con- 
currence matrices is up to ""("^H" -1 ) wn ich is greater 
than r g for all the systems greater than two qubits, a 
qubit and a qutrit and a qubit and a four-level system 
(a clear indication that the bigger the space, the more 
exceptional separable states are, as it was suggested by 
numerical analysis in j2||). To identify the variety of so- 
lutions of the equation set (|ll|l is the most costly part of 
the procedure. 

The necessary variable elimination involved in the pro- 
cess and based on a standard construction of a Groebner 
basis is known to be computationally difficult j2^| and 
even for 0-dimensional ideals, it may be polynomial in 
d n , where d is the maximum degree of the generators 
(original polynomials) and n is the number of variables 
[28|. This represents a super-exponential growth and, 
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like a brute force method indicated in the introduction, 
it may lead to prohibitive calculations already for small 
systems. 

Fortunately, by employing the XL (extended Lin- 
earization) method 29], one may drastically reduce the 
level of complexity with growing A = Af — r + 1 . The 
method relies on the observation that the set l|14|l can be 
regarded as a set of Af linear equations for ( r '£ 1 ) variables 
x a p — u a Ufj. Moreover, by multiplying each of the equa- 
tions by all the possible products u\ 1 • • • uj? r , such that 
Y^,iPi = D — 2, one can expand the original set into a ho- 
mogeneous set of Af( r ^>^~ 3 ) linear equations for ( r+ p~ 1 ) 
variables x 01) ,,, aj3 = u ai • • • u aD . The task is now to 
make the number of equations big enough to eliminate all 
but the last D + 1 variables, e.g., uf 3 , uf~ 1 U2, ■ ■ ■ , u§ 
and put, e.g., u\ = 1, so that one is left with a single 
polynomial equation in one variable. 

Numerical tests of the XL method reported in |29j in- 
dicate that for A = 0, the degree D of the single variable 
polynomial necessary to solve in order to obtain the solu- 
tion set of 1|14|) grows with r like 2 r_1 . This puts bipartite 
separability of generic states in 3 x 3 and even in 4 x 4 
dimensions within a reach of a simple PC. Actually, it 
may be even simpler than this. In 3 x 3 dimensions, a 
generic r is 9 just like a generic Af, thus giving A = 1. 
Numerical tests reported in [2j| indicate that in this case 
D = r which is clearly manageable even for relatively 
large r. Moreover, according to the same source, for 
A = 2 one has D y/r and for Af — er 2 , one should 
expect D « \\/y/e\. This indicates that for generic cases 
in many dimension (Af ~ r 2 /4) one should not expect 
the numeric difficulty of the problem to grow substan- 
tially with the dimensions of the local Hilbert spaces. 

When equation set iQ is solved then, to decide ex- 
istence of the required isometry, one should arrange the 
possible rows (ufei, Wfc2, "fc3, ■ ■ •) normalized by suitable 
factors Afe , so that the resulting matrix U is an isometry. 
For that, the factors have to satisfy: 

^\\k\ 2 u kt u* kj = (15) 

fc 

This is a simple set of linear equations for |Afc| 2 and it 
can be solved (separability) or proved unsolvable (entan- 
glement) without any problem. 

III. BEYOND BIPARTITE STATES. 

Our strategy easily generalizes to multipartite states 
with arbitrary number of parties sharing the state. One 
only has to analyze eigenvectors of the sum of all the in- 
dependent positive matrices JL l j£ defined in J2jj instead 
of the eigenvectors of a biconcurrence matrix B^ v dis- 
cussed earlier. In general, this will increase the number 
of concurrence matrices to consider, thus making full sep- 
arability less likely. Everything else goes like for bipar- 
tite states. For instance, for a 3-qubit state, there are 



three matrices A associated with the three types of pos- 
sible bi-partite correlations referred to in [2(j as , c| 
and 4 3) . The maximum possible rank of each of the 
matrices is three (the symmetric projector is three di- 
mensional). In a generic three-qubit state, one may then 
expect nine concurrence matrices (Af — 9) and r g = 8. 
To reduce the number of variables to one in the orig- 
inal equation set (|14|l one will then need such D that 
(t+d-ij _ D < A/"(''+ D 2 3 ). Here D = 5 satisfies the in- 
equality. In general, the bigger the system, the lower D 
will be sufficient. The difficult part of the algorithm is 
in the number of the (linear) equation and variables to 
consider. In the example above, one should expect them 
to be 1080 and 792 respectively. 



IV. CONCLUSIONS. 

We have pointed out that separability problem for ar- 
bitrary multi-partite states can be fully solved by a finite 
size, elementary recursive algorithm. However, in the 
worst case scenario, the underlying numerical procedure 
may grow super-exponentially with the state's rank. 

In an attempt to reduce this complexity, we investi- 
gated how analysis of concurrence matrices of given bi- 
partite and even multipartite states may lead to exhaus- 
tive analytic separability checks. It appeared that for 
generic states in m x n dimensions, analysis of concur- 
rence essentially reduced the task of solving separabil- 
ity problem to solving a set of linear equations in about 
/mn+c-ij variables, where D decreases with ran and for 
large ran it should not exceed 4. Moreover, the the same 
method is also applicable to multipartite states where it 
is at least equally efficient. 

One can notice, however, that the relatively low ex- 
pected complexity of the analysis of the generic states 
may be nothing more than a reflection of the fact that 
these states are usually entangled. Indeed, when the va- 
riety of the solutions of l|14H is either zero-dimensional 
or empty, then there is not much room for i|15fl to have 
any solutions. One possibility to extend the present re- 
sults can then aim at gaining some understanding of the 
intermediate case when there are many independent con- 
currence matrices, nevertheless the variety of solutions of 
(|14f> is at least one-dimensional. It would also be inter- 
esting to see whether in the language of concurrence it 
is possible to perform decomposition of the state into a 
separable part and a low dimensional part like it was for 
the edge states in optimization process ^3|. A possibil- 
ity like this would realize a version of Lewenstein-Sanpera 
decomposition |3Cj . 
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